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Abstract 

In this report, we present the derivation of the differential reflection spectrum as has been 
reported in Phys. Rev. B 72, 195301 (2005). 
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The differential reffection spectrum is derived based on a multiple layer structure with 
identical QDs, as schematically depicted in Fig. [TJ We will consider InAs QDs embedded 
in GaAs with an optical transition energy of the ground state far below the GaAs band gap 
energy (1.52 eV at 5 K). This assumption is well justified for InAs QDs, which have ground 
state transition energies of 1.1 eV or less The interaction of radiation and free carriers 
described by the Drude model, which is important at terahertz frequencies, is neglected in 
our derivations. 

For the derivation of the differential reflection equation, a procedure similar to that of 

n 

Tassone et al. [2] is used. Although, this method was initially developed to describe the 
response of quantum well systems, we will assume that the model is also suitable for the 
description of the reflection response of QD nanostructures. Justification of this assumption 
lays in the fact that the QD size is small compared to the wavelength of the probe laser. 
That is, the QDs can be treated as electrically small objects 0, 0]. The electromagnetic 
response of the nanostructure can be derived by means of an effective medium theory jj]. 
Because the QD volume fraction fqc is significantly smaller than the volume fraction of the 
host medium, we will use the Maxwell- Gar nett approach 0,0,1^1 to describe the dielectric 
properties of the nanocomposite. In this approach, the effective permittivity Seff of the QD 
layer is given by 

£e// - _ j,^^ £qd - £h 



with Eqd {sh) the permittivity of the QDs (GaAs host medium). For the calculations as 
presented below, a QD layer will be approximated by a thin layer with an average refractive 
index, uqd = y/e^JJ- 

If an electromagnetic wave propagates across an interface between two media with dif- 
ferent optical properties, its intensity will be divided between a refracted and a reflected 
wave. Assume that the wave propagates from vacuum into a dielectric medium with com- 
plex refractive index N = n + in, with n the real refractive index and k the attenuation 
index. The boundary conditions arising from Maxwell's equations require that the tangen- 
tial components of the electric and magnetic part are continuous across the boundary. Using 
the boundary conditions, the reflection of the electromagnetic wave can be presented by the 
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FIG. 1: Illustration of a QD nanostructure with n layers of QDs. The probe light is partly 
reflected and partly transmitted at the air-GaAs boundary. Subsequently, the transmitted light 
will be reflected by the different layers in the structure, and contribute to the total reflection. The 
order of the reflection is determined by the number of reflections within the structure. 
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cos9 - {N^ -sin^e) 
cosd + {N^ -sin^y/^' 



(2) 
(3) 



with r„ and Vp the reflection coefficient for the electromagnetic field components perpendic- 
ular and parallel to the plane of incidence, respectively, and 6 denotes the angle of incidence. 
The transmission coefficients are given by t„ = r„ + 1 and Ntp = + 1. The attenuation 
index of the dielectric medium causes a phase shift, and hence, r„ and Vp become complex. 
In case of normal incidence 6 = 0, the complex reflection coefficients are written by: 

n + in — 1 



n + in + 1 
3 
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with S the phase shift gained by the wave, which is given by tan(5= 2n/{n'^ + — 1). 
Note that also the transmitted wave has a phase shift. The intensity of the transmitted 
electromagnetic field is given by [l^ : 

with Jo = l-EoP the field intensity before penetration into the dielectric medium. R denotes 
the reflectivity given by R= |rp, a is the absorption coefficient of the medium, and d is the 
penetration depth. 

A schematic illustration of the QD nanostructure and the reflections induced by each layer 
is depicted in Fig. ^ The angle of the light path with respect to the QD plane is only by way 
of illustration. Hereinafter, the electromagnetic wave is assumed to propagate perpendicular 
to the QD plane. The total light reflection due to the QD structure is a superposition of the 
reflections of all interfaces. Taking into account all interference terms between the individual 
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ace reflections, the total sample reflection is calculated analytically in 
3], and can be described as 

\rs + CkTQDe'^' + C,ne'^^ + ^ D^r^rl^e^'*" + ■■■ + M\\ (6) 

k k 



Where denotes the surface reflectivity = "^''^^ , with ncaAs the refractive index 
of GaAs. rQ£) and represent the reflectivity of a single QD layer and the back surface, 
respectively. The summation takes into account the contributions of the individual QD 
layers with a relative phase 0fc induced by the optical path length of the GaAs barrier 
layers. C and D are correction terms which take into account the amplitude change due 
to reflection and absorption of the different layers, and M represents the higher order and 
mixed reflection terms. That are the terms due to multiple reflections on different QD 
layers. The first three terms on the righthand side of Eq. (jH)) are the zeroth order reflections 
(single reflection) whereas the double reflection on the QD layers is taken into account by 
the fourth term, etc. In this model, we will neglect the absorption of the probe light within 
the GaAs barriers, because the absorption coefficient a of GaAs is small within the energy 
window of the probe light , and because the GaAs layer thickness d is usually 

relatively small. Hereby, the amplitude change is primarily determined by the surface of the 
sample with transmittance tg. In addition, for GaAs the value of the attenuation index is 
much smaller than the real part of the refractive index [k <^ n) within the probe energy 



window. This means that, we can neglect the induced phase shift: S = 0. We can write for 
the correction terms Dc^C = t1 = l — r"^, for all k. 

Within the experiment, discrimination between the surface reflection, including the QDs, 
and the backside reflection can be made by changing the position of the probe laser focal 
point. Therefore, we will neglect the contribution of the backside to the total reflection 
signal. We can rewrite Eq. without loss of generality, as 

^ 1 - CrsrQDC'^" 
Due to the fact that rgD <^ 1, Eq. ((7j) can be simplified by 



k 

hence, 

k k k 

lin homogeneous QD layer in the vicinity 



For normal incidence, the reflection coefficient of a t 
of its transition energy (hujo) can be written as 
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TQDiuj) = ^ , ..-p , T, (10) 

where HT denotes the radiative broadening and hj the homogeneous nonradiative broadening 
of the reflectance spectrum. Rewriting Eq. (fTn|) as 

/ N (r + 7) + i(a; - Wo) p .^^s 

= - (^-^„)2 + (rH-oO' r. (11) 

and introducing the QD reflectivity into Eq. 0, we obtain 

+cV^.r^ ''"'^^° " \ (12) 

COS(0o) - 1 

Here, we assume z layers of QDs with equal spacing dgu, which introduces a phase difference 
00 due to the optical path length between adjacent layers. Hereby, the phase of each QD 
layer becomes (pk = k(f)Q. The third term in Eq. p2|) is small with respect to the first and 
second term, and will be neglected hereinafter. Equation (|12p can be simplified to: 

Rtot{io)=Rs{iu) + RQD{io). (13) 



The total reflection becomes a summation of the surface reflection Rg and the QD reflection 
Rqd- In this case the expression for the QD reflection can be written in the form Rq^iuj) 
= C{u!)T{u!), by introducing a line shape factor C{uj). Note that we postulate a frequency 
dependence of F. That is, F is assumed to be strongly dependent on the QD transition 
energy. 

In two-color pump-probe differential reflection spectroscopy, the photon energy of the 
excitation pulses is tuned above the GaAs band gap energy. Within the probe energy window 
which is far below the GaAs band gap, the pump induced changes of the GaAs refractive 
index can be neglected I3- Hereby, the refractive index changes are considered to be 
instantaneous with respect to the carrier occupation dynamics within the QDs. Therefore, 
the change of the sample dielectric function can be ascribed to the carrier induced change of 
the QD absorption, da/drj, as a result of energy level occupation in the quantum dot. Here, 
T] is defined as the carrier density. In addition, due to the QD energy level occupation the 
radiative broadening will be altered. Thus for two-color TRDR, the pump-induced reflection 
change within the probe energy window can be written as 

l^-(-) = ^Rqo(-) = • |.[^(-)r(-)], (14) 

which can be rewritten to: 

= -AT{lu)[T{lu)C'{lu) + C{u)], (15) 

with AF = ^ and C = -^jC. Note that the phase in the contribution of the QD layers is 
defined with respect to the sample surface reflection [Eq. (fT^ ]. but can now be interpreted 
as a relative phase between the individual layers. Moreover, in case of a symmetrical con- 
figuration, see Fig. [H the odd parts, i.e., the terms with sin0fc, will cancel. Hence, C{uj) 
becomes a typical Lorentzian line function. Although, we have made some assumptions, Eq. 
()15|) is a general expression for the carrier induced reflection change of a plane of identical 
QDs independent of the sample configuration. 

For realistic QD nanostructures, due to the QD-size distribution a summation over all 
homogeneously broadened QD states N{uji) has to be taken into account. This results in an 
inhomogeneously broadened QD density of states (DOS) of the ensemble, denoted by D{u), 



6 



such that 

Y^^RtoM) ^ - J][Ar(^,)H(^i)]iv(^i). (16) 

i i 

The terms in brackets can be interpreted as a response function of the QD reflection due to 
carrier occupation of the QDs, with Ti.{uj) = T{lj)C'{uj) + C{uj). Hence, the total reflection 
change becomes 

^Rtot{uj) = -n{Lu)D{iu), (17) 

with 

/oo 
[AT-n- N]duj. (18) 
-oo 

Due to the finite width of Ti, it is more appropriate to substitute ^ ATTiN by j[ATl-LN]. 
From Eq. (|T7jl we conclude that the differential reflection signal is determined by the 
QD DOS and the QD response function Til, which includes the carrier lifetime Td = F"^, and 
the carrier induced change of the radiative emission rate AF. 

This work is financially supported by the Dutch Foundation for Research on Matter 
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